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TECHNICAL MEMORANDUM 1425 

ON THE CALCULATION CF SHALLOW SHELLS* 
By S. A. Ambartsumyan 


1. We shall consider a sufficiently thin shallow shell of nonzero 
Gaussian curvature. In this case, neglecting certain small magnitudes, 
the problem, as shown by V. Z. Vlasov (ref. 1), can be reduced to a system 
of symmetrically constructed differential equations as follows : 

gg vfvf # - (HVi - LVh) V - 0 

- (HV§ - LVg) i - v|v| V + Z = 0 (1.1) 


These equations were constructed by the mixed method through the intro- 
duction of only two functions, namely, the stress function $, and the 
displacement function w. 


The forces T^, Tg, and S are expressed through <£ as follows: 

_ 1 3 /l S§\ 1 _ 1 d /l d§\ 1 9 a d# 

T 1 * B \B + T Z " A BE \A So/ + ^ 3p 3p 
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- Sr 


- — (■ 
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dcxdp ” A 3a ~ B 3p 3ay 


( 1 . 2 ) 


For the transverse forces N-j_ and Ng, we have the formulas 
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12(1 - v ) 
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(1.3) 
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In these formulas, 8 denotes the constant thickness of the shell, 
v the Poisson coefficient, E the modulus of elasticity, A m A (a, (3) and 
B » B(a,0) the coefficients of the first quadratic form of Gauss, 
ki =» k^(a,0) and kg =» kg(a,0) the principal curvatures of the coordinate 

surface in the orthogonal coordinate curves, 0 = constant, and 
a = constant. Further, the differential operators of ..the second order of 
the elliptical and hyperbolical type are defined as follows: 


. * (k b\1 

+ S 0 \B 


e AB I 3a \A Ba/ 

v h [ b2 h, (s“ sk) * a2 It (m ^)] 

The mixed operator HV| - LV 2 , in which H » l/2 (k-j_ + kg) and 
L ■ l/2 (k-j_ - kg), is defined by the formula 


(1.4) 


K - LV h) a ® [be (l k 2 ss) + ^ (i k i s^)] ( 1 - 5 ) 

2. We present equations (l.l) in a somewhat different form, so that 
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JL_ 

ES 




» vr * - vt w 




0, 


- V k $ - DV e w + Z » 0, 


( D 


E8 3 


12(1 - v 2 ) 


( 2 . 1 ) 


As shown by B. G. Galerkin (ref. 2), we note that the first equation 
of the system (2.1) is satisfied by introducing a certain displacement 
function cp (a , 0 ) , through which the required unknowns of the system (2.1) 
are expressed as ' > 




$ 


E8vJ<p 


(2.2) 


The second equation of equations (2.1), then assumes the form 

-_j2' 

6^ 




(2.3) 


Considering equation (2.2), we give the formulas expressing the 
computational magnitudes in terms of the displacement function as 
follows : 


E8 

T i - IT 


Tb / 1 1 sb s~] y2 E8 ra (1 s\ 1 ba b~| 

[Bp \b Spy BE Bo) V** t 2 - T [BE \a Bay + P' Bp Bpj 

E8 r_£_ 1 BB B 1 BA b! o2 

AB [dadp ” A Ba Bp B Bp BaJ k 


7 k<P 

(2.4) 
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fl d fl d\ 1 dA d v d /l v dB d”| o4m 

1 “ D |a 3a \A 3a/ + ab2 3p 3p + B 3p 3p/ + A 2 B cE e 

_ fl d fl d\ 1 dB d v d fl d\ v dA d""l V 4 W 

2 - D |b 3p \b 3p; + 3S 3S + A 3£ \k 3a/ + ^2 3p SpJ Ve<P 

_ _ D /- % r d 2 1 5 b S 1 dA d”! o4 

H i = - E 2 = - ab ( x - v ) [35Sp " a 3£ 3p " I 3p 3a| 


Nl 


_ D d „6 

a 3£ Ve<p ' 


N 2 


-2 & V°, 

B 3p e 


Thus, the problem of computing shallow thin shells with arbitrary 
normally applied loads reduces to finding the displacement function 
q> » > which is determined by the differential equation (2.3). 

3. Investigations by Y. N. Rabotnov (ref. 3) and A. L- Goldenveizer 
(ref. 4), show that the coefficients of the first quadratic form A and 
B, for a certain part of an arbitrary shell, behave almost like constants 
in differentiation. Hence, in the differentiation of products of the 
form Aw (or kw), the derivative of A may be neglected, and we can 
set d(Aw) * Adw. 

If it is assumed that a, and p are absolute coordinates, then, on 
' the basis of the previous discussion, we may set A * B * 1; we then 
obtain 


£ 

da 


d 2 


dp 


*2 


da £ 


+ k l (3.1) 


dp 2 


Equation (2.3) then takes the form 


_8 . 12(1 - v 2 ) n 4_ Z 

^ <P + ^5 V k ‘P ■ d 


(3.2) 


Now equations (2.4), (2.5), (2.6), and (fig. l) , are determined hy 
the formulas 


d 2 2 

T i - ES 'fr* 
^2 2 

T 2 = E8 7T v k^ 
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(3.3) 


“ ' V > 3S3p ^ ° 2 " D f^ + V |^ V4<P 
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'1 - - D & V 6 *, 91 " N 1 + § 


= - D 


[1; V2 


+ (1 - v) 


• (3.4) 


H 2 - - D |j v«*, % - »2 4 s - * D ' 2 + ^ - v > ^ 


4. As an example, let us consider a shallow shell, rectangular In 
the plane, freely supported on its contour, and subjected to a normal 
load. Let a and b denote the dimensions of the shell in the direc- 
tions a and p. _ 

The boundary conditions of the problem are 


w = 0, u = 0, T]_ = 0, Gq = 0 for a = 0, a = a 

w = 0, v = 0, T£ = 0, Gg = 0 for 3 = 0, 3 = b 


(4.1) 


These boundary conditions are satisfied by a solution of the form 


00 00 

\ ^ . . m ra . nrtB 

<P - / ,/ ,Ann s in — sin 


(4.2) 


Expanding the external normally applied load into a double trigonometric 
series, and substituting the results obtained in equation (3.2), yield 


4a° 

Am »~r 


a pb 


Z sin 2^2; S i n ^ dp 
a b 


0 UO 


(4.3) 


where m and n are odd positive numbers. 


Ain m [(m 2 + X 2 n 2 ) 4 + CxjJ (a > b, X - a/b ) (4.4) 


K = ^=^> C ~ 12 


(4.5) 


5. Let us consider the loading of the shell by a concentrated force 
P, applied at an arbitrary point (tj,£). According to equation (4.3), we 
have 


Am - - f a sin 223 sin 225 

TWrSo-hAI a * b 


(5.1) 
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Hence, from equations (4.2) and (3.3), we have 


<P = - 


T 1 - " 


4Pa^ 

Dit 8 ab 


m n 

4P\ 3 R 1 

a 2 2 — '2-—/ A mr 


(5.2) 
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\an 

mn 


^mx®np^mT)®n? 


^rax^np^mTi^n? 


m n 


S 


4PX 2 R 1 
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In the present and succeeding discussion, there is introduced for 
"brevity, the notations 


sin -f- =* Smx, sin =■ S^j 


mica 

a 


COS S 3Cro , COS 




sin 2^3. = Siqtjj 


cos = 0,^^; 


sin = Sn? 


COS 


rutC 


(5.4) 


= c n£ 


Further, from equation (2.2), we obtain 


w = - 


4Pa 

4 

Die ah 


YY ^ aE " 1 s, “ SnpSm ’' Sn 


(5.5) 


m n 


Making use of the identity 


EE 


(m 2 + \ 2 n 2 ) 


m n 


A« 

mn 



m n 


(m 2 + \ 2 n 2 ) 2 


EE 

m n 
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where 


(m 2 + X 2 n 2 ) 2 


(5.6) 


(5.7) 
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we obtain from equation (5.5), 


w 


4Pa ft 

Drt^ab 


551^ * Si " 554 -**■** 


The first term of this formula is the expression for the deflection 
of a rectangular plate, freely supported on its contour at the sides a 
and b under a concentrated load P, and applied at an arbitrary point 
Denoting the deflection of the plate by w*, we obtain 


w = V 4 cp 


* _ 4Pa* 
D* 4 ab 


C y y x»i — I (5.8) 

m" / '‘ TT / mn 1 


J 


The possibility of separating the deflection of a plate from the 
expression for the displacement w of a circular cylindrical shell was 
first shown by T. T. Khachatryan, in an unpublished doctoral dissertation. 

Further, making rose of the new expression (5.8), we obtain 

- 0.2 
> A -rnr 


°1 - 4 - C " 2 +vx2n2) 

°2 - 4 - jr yyl + vmZ> ‘www 5 * 

m n 

E - E * - (1 - U > 0 yy^ m <WWW>n£ 

m ti * nn 

N i “ ®i ~ Sr c ^y)y ^ m fr 2 + x2 * 2 ) c^ps^s^ 

N 2 = - Sr c^T> n (m2 + x2n2) WbpWbj 

' m n 11111 

%. - «f - fr °yy 'h: m + vxV > <WW 3 ®l 3 n« 

in n 

02 = §2 - ZfT' 0 ^ n ^ X2n2 + vm2 ^ ^CnpSmnSn^ 


(5.9) 
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In the preceding formulas, the magnitudes denoted with asterisks 
correspond to the tending of a rectangular plate (axb), coinciding with 
the contour of the shell. For computing these magnitudes, the well- 
known tables of B. G. Galerkin (ref. 5) may he used. 


The r emaining magnitudes in formulas (5.9) represent the effect of 
the additional internal forces that arise as a result of the curvature 
of the shell (the curvatures k]_ and k 2 enter in equation (2.1) in 

the composition of the secondary differential operator of the mixed type, 
and play the part of coefficient of elastic base, with respect to the 

flat plate). 

By an analogous method, the formulas may be obtained for a dense, 
uniformly distributed load q. = constant, over the entire surface of the 
shell, and are shown as follows: 


16 qa 8 \ ' s mx s n6 

D* 10 4-4— Ain 


Tl= - 


16ql 2 H 1 


c / / AdoAiP 

/ /Z— . — 

m n 


16qR ! \ 

Itr 0 2 /2 ^ 
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16 qXR]_ \ ^ \rm 
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% m n 
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6. In conclusion, a numerical example is given of the computation 
of a shell, rectangular in the plane, and freely supported over the 
contour for a uniformly distributed normal load. 

We take X = a/b = 2, the ratio K «= R;]/R 2 = 2, the Poisson coeffi- 
cient v a 0.3, the intensity of the load q, and Young's modulus E. 
According to equations (4.5), C = 20. 


For the deflection w, we obtain the formula 

.3 


w 


(w* - Aw) 


E5' 


( 6 . 1 ) 


where w* is a numerical coefficient taken from the previously mentioned 
tables (ref. 5), and Aw is an additional deflection depending on the 
curvature of the shell, so that • 


Aw a 


16x12(1 - v 2 ) C\ 4 ^mn mica. 




sin — — sin ~i. 
a d 


(6.2) 


Restricting ourselves to two terms of the expansion in formula 
(6.2), the following results of the computations for several points (<x,p) 
of the shell are given: 


Points 

(a - i a, p - \ ^ 

^a, a 

\ a , 0 - \ b ) ( a 

*“ g" a , 3 ™ 


w* 

0.0843 


0. 1106 

0.0760 


Aw 

.0449 


.0608 

.0431 


w* - Aw 

.0394 


.0498 

.0329 


For 

the bending moments 

Gl 

and G2) we obtain the formulas 


Gi « qb 2 (m* - 


, G 2 = qb 2 (m£ - 

An 2 ) 

(6.3) 


where m* and m^ are coefficients taken from table III (ref. 5), and 


16\ 2 C \ N \~ A lmn(m 2 + vX 2 n 2 ) . mjta . rot [3 

^ — S3!: sln — Bln T 6 

m n 


Ang 


16X 2 C \ N \ v Xmn(X 2 n 2 + vm 2 ) 

sin mjta 




m 


TrmA F 
mn 


sin 222 sin 

a b 


(6.4) 
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Restricting ourselves to two terms of the expansion in formula 
(6.4) , the following results of the computation for several points are 
given : 


Points ^c. = 

T a, P - | b) (a - J 

a, («• 



■? 

0.0446 

0.0464 

0.0304 


An! 

.0241 

.0276 

.0195 


m£ - An^ 

.0205 

.0188 

.0109 



.0807 

.1017 

.0773 


Ai^ 

.0444 

.0583 

.0412 


m| - Aug 

.0363 

.0434 

.0361 


For the forces 

T! and Tg, we have 

the formulas 




Ti 53 * <3?iAt!, 

Tg = - qRlAtg 


(6.5) 


where 


At x 



. mita . 

sin sin 

a 


nitp 

h 


At 2 


16C 

1? 



. mm . 

sin sin 

a 



Restricting ourselves to three terms of the expansion in these 
formulas } the results of the computations are: 


§*) f 


At^ 

0.4252 

0.5569 

At 2 

.1333 

.1010 


Points 


(& - J P " \ (a- 



0.3928 

.0711 
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